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(N : ABSTRACT 



During a Hubble time, cluster galaxies may undergo several mutual encounters 
Q> i close enough to gravitationally perturb their hot. X-ray emitting gas flows. We ran 

several 2D, time dependent hydrodynamical models to investigate the effects of such 
perturbations on the gas flow inside elliptical galaxies, focusing on the expected X-ray 
properties. In particular, we studied in detail the modifications occurring in the 



0^ 



Ph' 



^ ■ scenario proposed by D'Ercole et al. (1989), in which the galactic interstellar medium 

c/3 ■ produced by the aging galactic stellar population, is heated by type la supernovae 

(SNIa) at a decreasing rate. We find that, although the tidal interaction in our models 
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lasts less than 1 Gyr, its effect extends over several Gyrs. The tidally induced turbulent 
flows create dense filaments which cool quickly and accrete onto the galactic center, 
C^ I producing large spikes in the global X-ray luminosity Lx- Once this mechanism starts, 

it is fed by gravity and amplified by SNIa. This evolution is found to be virtually 
independent of the dynamical state of the gas flow at the beginning of the interaction. 
To better understand the role of SNIa heating, we also considered a "pure" cooling 
flow model without supernovae; in this case the amplitude of the Lx fluctuations due 
to the tidal interaction is substantially reduced. We conclude that, if SNIa significantly 
contribute to the energetics of the gas flows in ellipticals, then the observed spread in 
the Lx — Lb diagram at any fixed optical galaxy luminosity Lb ^ 3 x W^^Lq may be 
caused, at least in part, by this mechanism. On the contrary, tidal interactions cannot 
be responsible for the observed spread if the pure cooling flow scenario applies. 



Subject headings: Galaxies: cooling flows - Galaxies: elliptical and lenticular, cD 
Galaxies: interactions - Galaxies: ISM - X-rays: galaxies 



1. INTRODUCTION 

In a previous paper D'Ercole et al. (1989) proposed the wind-outflow-inflow (WOI) scenario 
as a possible explanation for one of the most striking properties of the X-ray emission of elliptical 
galaxies, i.e., the large scatter in the Lx — L-q diagram of roughly two orders of magnitude in 
Lx at any fixed Lb ^ 3 x 1{)^^Lq (Fabbiano 1989, Fabbiano, Kim & Trinchieri 1992). In their 
subsequent extensive exploration of ID hydrodynamical models Ciotti et al. (1991, hereafter 
CDPR) assumed that the type I supernova (SNIa) rate decreases in time as i?sN oc t~^'^, faster 
than the decrease of the mass return rate from the (passively) evolving stellar population, which 
varies in time approximately as t^^'^^. In these models, the gas lost by the stars is initially driven 
out of the galaxy through a supersonic wind powered by the thermalization of the kinetic energy 
injected in the interstellar medium (ISM) by the SNIas. As the specific energy input decreases, the 
wind turns into a subsonic outflow and eventually the flow reverts to an inflow regime, at the time 
of the so-called cooling catastrophe (tec)- In pace with the content of the X-ray emitting hot gas, 
Lx decreases during the wind phase and increases during the outflow regime up to a maximum 
reached at tec- Later on, Lx slowly decreases roughly following the quenching of the SNIa heating 
LsN oc i?sN- CDPR ran several ID numerical simulations, showing how tec critically depends on 
the values of the structural parameters of the host galaxy, such as the scale-length of the stellar 
distribution (e.g, the effective radius for de Vaucouleurs proflles or the so-called core radius for 
King profiles), the scale-length of the dark matter halo, the relative amount of total masses of the 
two distributions. Small changes in one (or more) of these parameters produce a large shift (of 
the order of some Gyr) of tco causing two galaxies of virtually the same Lb to be caught at very 
different Lx after a Hubble time. As a consequence, the above-mentioned Lx scatter at fixed Lb 
finds a simple explanation in the WOI scenario. 

The main objection raised against this scenario has been the claimed detection of low iron 
abundance in the ISM of ellipticals. In fact, under the assumption of solar abundance ratios, 
the analysis of the available data suggests a very low iron abundance, consistent with no SNIas 
enrichment, and even lower than that of the stellar component (Ohashi et al. 1990, Awaki et 
al. 1991, Ikebe et al. 1992, Serlemitsos et al. 1993, Loewenstein et al. 1994, Awaki et al. 1994, 
Arimoto et al. 1997, Matsumoto et al. 1997). However, some authors have found that more 
complex multitemperature models with a higher abundance give a better fit to the data (Kim et 
al. 1996, Buote & Fabian 1998). This is rather puzzling, in that a value as high as one fourth 
of the standard Tammann (1982) rate i?sN = 0.22 SNU (for Hq = 50 km s^"*^ Mpc^^) agrees 
with the current optical estimates of the present day SNIa rate (Cappellaro et al. 1997). Several 
models have been calculated with this new rate (Pellegrini & Ciotti 1998, D'Ercole & Ciotti 1998, 
hereafter DC), resulting in the main features of the WOI scenario being preserved, provided the 
amount of dark matter is properly scaled down. 

In any case, even though i?sN ^ 0, it is well known that the pure cooling fiow scenario is 
unable, in isolated galaxies, to reproduce the observed scatter (see, e.g., CDPR). As a consequence, 
the possible environmental origin of the Lx — Lb scatter has been proposed. Although some 



discussion was dedicated to environmental effects, tlie bulk of the CDPR simulations did not 
consider ambient influences. It is known, however, that the majority of early type galaxies is 
found in clusters, where they certainly suffer several environmental interactions. More specifically, 
1) thermal conduction with the intracluster medium (ICM); 2) ram pressure stripping of the ISM 
due to the motion of the galaxies with respect to the ICM; 3) influence of the cluster tidal field; 4) 
galaxy-galaxy gravitational encounters. 

Because the ICM in rich clusters is much hotter (T ~ 10^ K) than the galactic gas flows 
(T ~ 10^ K), heat is expected to flow into the galactic ISM if thermal conduction is not suppressed 
by magnetic fields. It is easy to show that the heat transferred by an unimpeded thermal 
conduction is much larger than any other galactic energy budget, and evaporation can hardly be 
avoided. The very existence of bright X~ray ellipticals thus suggests that heat conduction must 
be suppressed to a great extent (see, e.g., CDPR, and references therein). 

The ram pressure effect of the ICM on the ISM of the elliptical galaxies was discussed by 
D'Ercole et al. (1993) in the WOI scenario, and, e.g., by Balsara, Livio & O'Dea (1994, and 
references therein) in the cooling flow scenario. While numerical models show that ram pressure 
is very effective in stripping galaxies of gas, observations reveal that numerous cluster galaxies 
are strong X-ray emitters, thus showing their ability to retain their hot ISM. For example, Dow 
& White (1995), found that ellipticals in Coma are as X-ray bright as in Virgo. This situation 
appears rather puzzling, and certainly not satisfactorily understood. 

The mean cluster gravitational field acts on the structure of each galaxy and its gaseous halo, 
as a weak tidal field. A''-body numerical simulations (e.g., Ciotti & Dutta, 1994) show that its 
effect on the galactic stellar component is very small, consisting essentially in a re-orienting of the 
galactic major axis with respect to the cluster center. 

In summary, while some effort was spent by several authors in considering the influence of 
the ICM on the ISM of ellipticals, up to now no attempt was made to study the response of the 
X-ray emitting ISM of cluster ellipticals to the tidal interaction due to the cluster field and to a 
close encounter with a nearby galaxy. As we will see, in a typical cluster, encounters close enough 
to have non negligible dynamical effects on the ISM are likely to occur. 

In this paper we show the results of 2D numerical hydro dynamical models focusing on the 
influence both of the cluster tidal field and of a transient external gravitational field on the 
evolution of the ISM of ellipticals in the WOI and cooling flow scenarios. In §2 we describe the 
galaxy model, the hydro dynamical equations, and the assumptions adopted in our simulations. 
The results are presented in §3. A discussion and the conclusions are finally given in §4. 
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2. The Galaxy Model 

In order to better compare the present results with those obtained by CDPR, we adopt for 
the model galaxy the same bright and dark matter (DM) distributions considered by these authors 
for their basic model (the so-called King Reference Model, KRM). The stellar density distribution 
is a truncated King (1972) model 

and the DM density profile is a truncated quasi-isothermal model 

PhO 

1 + r"^ /rl^ 

Here (p*0)'''*) a-iid (/9h0)'''h) are the central density and scale-length of the luminous and DM 
density distributions, respectively, while r^ is their common truncation radius. The total galaxy 
density, potential, and mass are Pg = p* + Ph-, (kg = (k* + 4>h-, and Mg = M* + Mh. 

As for the KRM in CDPR we assume p^,o = 6.08 x 10^^^ g cm^^, pho = 7.6 x 10~^^ g cm~^, 
r* = 368 pc, Th = 4.5r* and r^ = 66.24 kpc. With these parameters M* = 2.75 x 10"*^^ Mq and 
Mh = 2.48 X 10^^ Mq, respectively, and the galaxy lies on the so-called Fundamental Plane. 
It should be noted that several recent indications, from high-resolution A^-body simulations 
(Dubinski & Carlberg 1991; Navarro, Frenk & White 1997) as well as from theoretical arguments 
(Evans & Collett 1997; Ciotti 1996,1999), seem to suggest a peaked density distribution for the 
DM halos, as for example that used in the numerical simulations of gas flows carried out by 
Pellegrini & Ciotti (1998). In spite of this, we decided to maintain here the functional form for 
the halo as in the KRM, in order to better understand the effects of the tidal interactions by 
comparison with the unperturbed model. 

Of course, a self-consistent treatment of the evolution of galactic gas flows in the cluster 
environment is a formidable task, both for the complexity of the involved physics, and for the 
high dimensionality of the parameter space (requiring very large computational resources to be 
fully explored). Our aims are far more narrow, focusing on the qualitative understanding of the 
effects of a representative tidal encounter. In this line, the mass Mp of the perturbing galaxy (the 
perturber) is assumed to be equal to the mass of the galaxy, Mp = Mg ~ 2.8 x 10^^ Mq, because 
we are not interested here in the effects of interactions with significantly less massive galaxies. In 
a subsequent paper (Pellegrini, D'Ercole & Ciotti, in preparation) we explore the effects produced 
on the X-ray emission of a galaxy model similar to that here discussed by a collision with a 
significantly less massive cluster galaxy. 



2.1. Assumptions and Equations 

We derive and discuss here the hydrodynamical equations describing the evolution of the gas 
flows in our models, using the general equations given in Appendix A. 2. In the following, bold-face 



symbols represent vectors. In the (baricentric) cluster inertial reference frame S, the center of 
mass of the galaxy and of the perturber are R and Rp, respectively; moreover, rp = Rp — R and 
Vp = llfpll, where || • || is the standard norm. For computational reasons these equations are written 
in a non-inertial reference system S' fixed by the following rules: the origin of S' coincides with 
R, ^ is the position vector in S', the ^3 axis points at all times towards Rp, and finally r = |||||. 



2.1.1. The Effective Gravitational Field 

We derive first the effective gravitational field experienced by the gas in the frame S' 
(see equation [A16]), and we discuss two assumptions used in our models, namely that 1) the 
cluster tidal field is negligible during the encounter, and that 2) the perturber trajectory can be 
approximated by a straight line, covered with constant velocity. Although in the numerical code 
the gravitational forces are not approximated by their linear expansion (tidal regime), here we 
consider the associated tidal fields, in order to obtain the geometrical conditions for the validity of 
the assumptions above. 

At any position x in the cluster frame 5", the total acceleration is g = — Vx(i;^ci + <Pp + (f'g)i 
where (pci, (pp, and (pg are the cluster, perturber, and galaxy potentials, respectively, and 
Vx is the gradient operator with respect to x. Note that the self-gravity of the gaseous 
component is not considered. The acceleration A of the center of mass of the galaxy is given by 
MgA = — /pgVx(0cl + (pp)^^^- Ii^ the tidal approximation cpci and (pp are expanded up to the 
second order around R, i.e., 

, / N , /^^ /^x ^ <rci(R)(x-R),x-R> 
0ci(x) ~ (/>cl(R)- < gcl(R), X - R > ^ '^ ^ '' , (3) 

where gci(R) = — Vr(/>ci(R), <,> is the standard inner product, and T^ is the tidal tensor 
associated with (pd, whose elements are given by 

A similar relation holds for the gravitational field produced by the perturber, where now 
gp(R) = — Vr0p(R) and Tp is the tidal tensor associated with (pp. In this approximation, simple 
calculations show that the effective gravitational field entering equation (A16) is given by 

g-A~[Tei(R) + rp(R)](x-R)+gg, (5) 

where gg = — Vx'/'g. The coordinate transformation from S to S' shows that 



0'gg = -G^^i, (6) 

and 






OT[rp(R)(x-R)] = -^(6,6,-26). (7) 



Equations (6) and (7) are derived under the assumption of spherical mass distributions. Actuahy, 
the galaxy and the perturber will change their shape during the interaction. However, for small 
deformations the potential remains rounder than the underlying density distribution (see, e.g., 
Binney & Tremaine 1987, hereafter BT), and so we regard equation (6) as still largely valid for 
our problem. Moreover, the effect of the deformation of the pertuber density influences the tidal 
tensor Tp given in equation (7) only at higher order in rp (see, e.g., Danby 1962, Chapter 5). As 
a consequence, in all simulations the (visible and dark) mass densities of the galaxy and of the 
pertuber are assumed to maintain their initial spherical symmetry. 

As anticipated in point 1) above, in our simulations we do not consider the term Td appearing 
in equation (5), in order to reduce the dimensionality of the parameter space. The influence of the 
cluster tidal field in absence of close encounters was explored separately, and it is discussed in §4. 
We now determine when the perturber tidal field Tp is dominant over Tc\. The cluster tidal field 
at distance R= ||R|| from the cluster center is given by equation (4), and its explicit expression is 



T^m = - T ' 10, (8) 




where 'Pd[R) = 3Mci(i?)/47ri?^ is the cluster mean density inside -R, and < q{R) = 
Pc\{R) /]m{R) ^ 1- The matrix expression of Td is obtained in the reference frame S with its z 
axis aligned along R (see, e.g., Ciotti & Dutta 1994). As a consequence, in equation (8) the radial 
and the tangential components of the cluster tidal field are apparentjj. In order to obtain a simple 
estimate of the relative importance of the cluster and pertuber tidal fields, we consider a typical 
cluster with a density profile pci as that given in equation (1), with a core radius r^ = 350 kpc, 
and a central velocity dispersion o"ci = 10^ km s~^ (Sarazin 1986). We compare the minimum 
(absolute) value of the perturber tidal field components, GMp/r^, with the radial and tangential 
components of the cluster tidal field. In Fig. 1 these quantities are showed, and it turns out that 
the cluster tidal field can be neglected for rp ~ 100 kpc everywhere but inside the cluster core. 



2.1.2. The Pertuber Trajectory 

Following the arguments above, we assumed in our simulations r™™ = 100 kpc as the 
minimum distance during the encounter, in order to safely neglect the cluster tidal force on 
the ISM. Actually, the (full) perturber gravitational field is taken into account for distances 
much larger than rjf^, i.e., for rp < 450 kpc. This in order to avoid a sudden change in the 
total gravitational field experienced by the gas flow, that could be at the origin of unphysical 
oscillations. In other words, up to the time t' (a free parameter in the simulations, whose choice 



^Note that the tidal field outside a spherical mass distribution given in equation (7) can be obtained directly from 
equation (8) for g = 0, and by substitution of pci with pp. 



is given in §3) the gas flow evolves as in CDPR, successively the perturber is placed at a distance 
of 450 kpc, with an initial velocity oi V = 1000 km s^^ relative to the galaxy. We indicate with 
tjn the time at which rp{t^) = r™™, and with Atp = tm ~ ^p- 

In order to follow the time evolution of the ISM, we obviously need to know how rp changes 
with time. Adopting equation (3) to describe the effect of the cluster gravitational field on the 
galaxy and pertuber motions, one obtains 

i^p ~ -G(Mg + Mp)^ + rei(R)rp. (9) 

With the assumed parameters for the galaxy, the perturber, and the cluster density distribution, 
it is easy to verify from the equation above that the effects of Td on the motion of the two galaxies 
are non negligible for values of rp in the range of interest, and so the relative orbit cannot be 
described by a two-body problem. As anticipated by point 2) in §2.1.1, we arbitrarily adopt a 
linear trajectory for rp, described at constant velocity V. In this assumption Atp ~ 460 Myr, and 
simple geometrical arguments show that the only non-zero component of the angular velocity in 
S", entering equation (A16), is given by: 
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We point out that, as a consequence of the geometry of the problem, the non-inertial forces 
associated with the rotation of S' cannot be taken into consideration by our 2D code, and so we 
drop these terms in the hydrodynamical equations used in the numerical simulations (equations 
[12]-[14]). In §4 we will calculate a posteriori, from the computed models, the ratio between the 
gravity and the fictitious forces associated with f2, in order to obtain a rough estimate of the 
severity of this approximation. Although this assumption may appear rather rough, the computed 
models show that the perturbations in the ISM do not depend on the details of the physical 
process which generate them. Thus our assumption is not particularly severe. 



2.1.3. The Source Terms 

As already discussed in the literature (Mathews & Baker 1971), the mass and energy 
sources in elliptical galaxies are associated with their evolving stellar population: the two main 
contributions are due to the gas lost by the red giants and to the ejecta of SNIas. The total mass 
return rate, expressed in the general case in equations (A15)-(A17), is given here by A^ = a(t)p*, 
where a{t) = a^:{t) + as^{t). We adopt the same a{t) as in CDPR (equations [3] and [7] there). In 
particular, this corresponds to a present day i?gN = 0.15 SNU, i.e., 2/3 of the standard Tammann's 
rate. Moreover, we assume that the internal energy associated with the mass ejection in both 
processes is negligible with respect to the thermalization of the stellar velocity dispersion and to 
the SNIa ejection velocity. The exact expression for the injection energy is given by 

„ a*u^ + asN^N , Tr(o-2) osn^n + «*Tr((T2) 
" 2a 2 2a ' ^ ' 



where the last formula is derived under the hypothesis that 1) the velocity of stellar winds is 
considerably lower than the stellar velocity dispersion (i.e., u^ « a"^), and that 2) the stellar 
velocity dispersion is lower than the velocity of SNIa ejecta (i.e., o"^ << u^^). In accordance with 
the previous works we adopt the second, approximate expression for £q, but we remind that in 
some circumstances the full one should be used. The cooling rate per unit volume is given by 
C = ncnpA(T), where we parametrize A(T) following the prescription of Mathews & Bregman 
(1978), after correction of a wrong sign in their equation (A2). Finally, no streaming velocity of the 
stellar component is assumed to be present in S", and accordingly in the momentum and energy 
equations (A16)-(A17) we set v' = 0. This assumption can be qualitatively explained considering 
that, when the mutual interaction between the galaxy and the perturber starts, their stellar and 
dark matter distributions are distorted. These deformations are responsible for the appearance 
of a net torque acting on the galaxy density distribution, from which an angular momentum is 
originated. Of course, the real situation is much more complex, but our simple assumptions should 
capture the kinematical behavior of the galaxy stellar and dark matter components. 



2.1.4- The Hydrodynamical Equations 

As a consequence of all the previous assumptions and discussions, equations (A15)-(A17) are 
specialized in the numerical code as follows: 

dp 

-Q^+y^- pi> = ap^, (12) 

dpu 



dt 



+ V| • {pu ®u) = p(gg + T^i) - (7 - l)V^E, (13) 

^ + V^.Eu = -{-f- 1)EV^ .u+^\\u\f + ap,£o - C. (14) 

p, E, and v> are respectively the density, the internal energy per unit volume, and the velocity 
of the gas, expressed in S'. The explicit expressions of gg and Tp| in this frame are given in 
equations (6) and (7). The gas pressure is p = (7 — 1)E, where 7 = 5/3 is the ratio of the specific 
heats. To integrate the set of equations we used a second-order, upwind, 2D numerical code 
described in DC, coupled with a staggered, spherical, Eulerian grid. Given the symmetry of the 
problem, the grid covers polar angles 0° < ^ < 180° with 82 equally spaced angular zones. The 
radial coordinate is divided into 72 zones and extends up to 95 kpc. In order to get a reasonable 
spatial resolution in the inner region, a geometrical progression is adopted for the radial mesh 
size, the first one being 100 pc wide. Reflecting boundary are assumed everywhere but at the 
outer boundary, where outflow conditions are adopted. With the adopted non-inertial reference 
frame 5" and numerical grid, the perturber is always outside the grid and moves along the ^ = 
direction (i.e., along the ^3 axis). As in CDPR, the gas temperature is not allowed to drop below 
10^ K, to avoid an excessive reduction of the time step when a rapid cooling is present. We assume 



the model galaxy to be initially devoid of gas due to the previous activity of the Type II SNe. A 
discussion on the reliability and implications of this assumption is given in CDPR. 



3. The Models 

Given the Courant condition on the angular mesh close to the center, the time steps for the 
numerical computations are rather short (few 10^ yr). We thus computed only few models. We 
first ran a model without the perturber and cluster gravitational field (hereafter model MO), i.e., a 
model analogous to the KRM presented in CDPR. As for the KRM, the cooling catastrophe occurs 
before 15 Gyr, and all the three dynamical phases - wind, outflow and inflow - are recoveredpl. 

In order to understand whether there are differences in the flow evolution depending on the 
time tm at which the maximum approach happens, we ran three models (Ml, M2, M3) for three 
different values of tm, namely t^^'^ = (6.3, 10.8, 12.6) Gyr. In each of these models the initial 
conditions are given by the hydrodynamical quantities of model MO at ti^ = t^ — Aip. In model 
Ml the galaxy starts to be perturbed when it is still in its wind phase, while in models M2 and 
M3 the galaxy is perturbed when is in the outflow and inflow phase, respectively (see Fig. 3). We 
also ran a pure cooling flow model (hereafter model MC) in which the energy injection by SNIa is 
absent. For this model we set f = 4.9 Gyr. Although this model is expected to be rather similar 
to model M3, we ran it to ascertain whether the presence of the SNe energy input may affect the 
perturbed gas flow, even in the inflow phase. 

Finally, we ran two more models without perturber, in order to check the sensitivity of our 
results to some of the assumptions made. Namely, 1) model MP is analogous to model M3, but 
with a random large-scale velocity field perturbation added to the initial conditions, as done by 
Kritsuk, Bohringer &: Miiller (1998) in their simulations; and 2) model MT is analogous to model 
MC, but with the cluster tidal field taken into account during a galaxy orbit inside the cluster. 



3.1. Models Discussion 

We discuss here in detail the Ml model, in which the galaxy starts to be perturbed during 
the wind phase. The other M2 and M3 models are presented by comparison with the Ml model, 
pointing out differences and analogies. 



^Note however that in model MO tec is greater than in the KRM (11 Gyr instead of 9.3 Gyr, see Fig. 3). This 
depends on the different value assumed for the mean molecular weight (0.62 instead of 0.5) which enters quadratically 
in the cooling term. By the way, this shows once more the extreme dependence of the WOI models on the parameter 
values. 
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3.1.1. Hydrodynamics 

Figure 2 shows the density contours and the velocity field of the gas at four different times. 
Panel a is a shot taken at a time very close to tm, when the distance to the perturber is slightly 
larger than r™™; the exact time and relative distance are (t,rp)=(6.4 Gyr, 141.4 kpc). The ISM in 
the galactic region facing the perturber accelerates toward the right direction (where the perturber 
is located), and is lost by the galaxy at higher velocity than in the opposite direction. This 
apparent asymmetry is due to the fact that in the code the full gravitational field of the perturber 
is considered, which, for short distances, differs significantly from the associated (symmetric) tidal 
field. In panel h [t = 9.9 Gyr) the effects of the encounter are still clearly visible, although the 
perturber distance from the galaxy (for the assumed rectilinear trajectory) is r^ = 3.6 Mpc. This 
is also true at later times (panels c and d), where eddies are established and different dynamical 
regimes are present, some regions being in outfiow and others in infiow, a situation analogous of 
that found in DC on their investigation on gas fiows in (isolated) SO galaxies. 



3.1.2. X-ray Luminosity and Gas Mass Evolution 

The X-ray luminosity evolution of model Ml is shown as a solid line in Fig. 3a. For 
comparison, the dashed line represents the Lx temporal evolution of the model MO, and the dotted 
line shows the time evolution of the SNIa energy input for unit time, LgN. 

As the encounter with the perturber occurs, the cooling catastrophe is anticipated, as a 
consequence of the increase in the mean density in the galactic central regions produced by the 
tidal perturbation, and Lx reaches its (first) maximum at tec — 9.4 Gyr, earlier than for model MO 
(ice — 11 Gyr). This is an interesting result. In fact, note that considering spherical coordinates 
with the ^3 polar axis (0 = 0°) oriented towards the perturber (the coordinate system actually 
adopted in the numerical code), from equation (7) one obtains that the radial component of the 
tidal force (for unit mass) is given by Frad = — (GMpr/rp)(sin^ 6* — 2cos^6') and is expansive for 
Q° < 9 < arctan\/2 ~ 54° and 126° < 9 < 180°. On the complementary spherical sector the 
field is compressive. Due to the strong non-linear dependence of the cooling function on the gas 
density, it is not easy to determine without the aid of numerical simulations which effect prevails 
in determining the successive evolution of the fiow, i.e., if the cooling catastrophe is anticipated or 
retarded by the effect of the tidal field. Our results show that actually compression wins the game 
and cooling increases. 

Successively, strong and fast oscillations in Lx occur, whose amplitude may reach an order 
of magnitude. These oscillations are due to correspondent oscillations of the ISM density in the 
galactic central regions: in fact, most of the X-ray luminosity is emitted there^, for Lx oc n^. It is 



^Being Lx ~ J n'^A{T)dV, the oscillations in Lx are basically due to variations and/or rearrangements in the hot 
gas density, because A(r) is nearly constant in the range of X-ray temperatures. 
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important to note that, while these density oscillations reflect dramatically on Lx, they have little 
effect on the overall content of the hot gas which decreases monotonically in this phase, as shown 
in Fig. 3d (medium solid line). 

This strongly unsteady evolution is due to the non-spherical gas accretion onto the galaxy 
central region. Converging to the center, non-radial gas streams interact and compress each other, 
giving rise to tongues of dense cold gas (see Fig. 4, where enlargements of the gas distribution 
and of the velocity field in the central region are shown) . These tongues are rapidly accreted onto 
the galaxy center. When the compressed gas quickly cools, Lx suddenly increases. The cooled gas 
disappears in the center, leaving behind regions of low gas density which do not radiate efficiently 
and tend to expand as they are heated by SNIa: as a consequence Lx abruptly decreases. New 
streams of hot gas coming from the outer regions contrast the expansion of the gas close to the 
center and form new cold filaments. This cycle is fed by gravity and amplified by Lsx. Actually, 
the hot gas mass decreases after every spike of Lx (Fig. 3d, medium solid line), while the cold gas 
mass increases specularly (Fig. 3d, light solid line). When the mass of the hot ISM is sufficiently 
reduced, Lx falls by nearly two orders of magnitude. Successively, the low density, non radiating 
hot gas is able to revert the fiow on large scale, and a new wind phase starts (see Fig. 3a and 
Fig.3d for t > 12 Gyr). 

The importance of SNIa in maintaining this cycle is apparent when considering the X-ray 
luminosity evolution of model MC, shown in Fig. 5. In this model, where Lsn = 0, the oscillations 
of Lx have smaller amplitudes and the dramatic drop shown by the model Ml is absent. In fact, 
the Lx oscillations in model MC never exceed a factor of three relatively to the unperturbed 
cooling flow solution. Concerning the total mass budget (Fig. 3d, heavy solid line), it is apparent 
its substantial reduction when compared to model MO (Fig. 3d, heavy dashed line). Thus, the tidal 
encounter is able to strip a significant fraction of ISM from the Ml model. 

Figure 3bc show that models M2 and M3 (where the maximum approach happens in the 
outflow and inflow phase, respectively) exhibit Lx oscillations similar to model Ml. It is however 
apparent that in these latter models the dramatic drop in Lx (occurring in model Ml at t ~ 12 
Gyr), is absent over an Hubble time. This is due to the fact that Lsn at the time of the encounter 
is lower than in model Ml (Fig.3abc, dotted line). Thus less energy is available to sustain a 
turbulent flow, and the Lx (and mass budget) evolution of models M2 and M3 is more similar to 
model MC. 

The hot and cold ISM masses of models M2 and M3 (medium and light solid lines in Fig. 
3ef) evolve in the same qualitative way as in model Ml (Fig. 3d). On the contrary, their total 
mass (heavy solid line) increases, at odd with model Ml. This result at flrst could seem surprising 
because model Ml is in a wind phase for t > 12 Gyr, and one could expect an easier degassing 
than in models M2 and M3. Actually, the hot gas density in model Ml is lower than in models 
M2 and M3 (as can be seen considering the amount of hot gas mass in the three models): as a 
consequence, the mass loss at the outskirts (r = rt) of the galaxies increases from model Ml to 
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model M3. In model Ml the mass loss rate is lower than the mass return rate from the evolving 
stellar population, a{t)M^:, the opposite is true for models M2, and this explains the different time 
behavior of the total mass in the three models. 



3.1.3. Surface Brightness 

The surface brightness profiles of all computed models do not result strongly influenced by the 
tidal interaction. As a representative case we describe the X-ray surface brightness distribution 
Sx of model Ml only. During all the model evolution, the isophotes remain rather circular, 
although locally distorted, even for a viewing angle 9 = 90° (i.e. a direction perpendicular to the 
line joining the galaxy and the perturber, where the effect is maximum). This can be seen in Fig. 

6, where two sets of isophotes are superimposed, referring to t = 9.9 Gyr (solid lines) and t = 11.1 
Gyr (dotted lines). These times correspond to a local maximum and minimum of Lx, respectively 
(black dots in Fig. 3a). Radial cross sections of Sx are shown (after an angular mean) in Fig. 

7. Note that when Lx is higher {t = 9.9 Gyr) the extra-luminosity is radiated mostly from the 
center, and so density oscillations here are by far more important than in the rest of the galaxy. 



4. Discussions and Conclusions 

In this paper we have investigated the consequences of a gravitational encounter on the 
dynamics of the hot. X-ray emitting ISM of cluster elliptical galaxies. We adopted a spherical 
model galaxy identical to the King Reference Model described in CDPR in order to better compare 
our results with the "standard" WOI picture. It turns out that once the gas is perturbed by the 
gravity of a galaxy passing nearby, it remains turbulent for the rest of the time. In fact, when 
the spherical symmetry is broken by the tidal interaction, non-radial streams converging toward 
the center merge thus compressing the gas, which quickly cools and "disappears" into the center. 
Some regions close to the center are then formed where the gas is rarefied and tends to expand as 
heated by the sourcess terms. This expansion is contrasted by the outer ISM which is falling to 
smaller radii, and new cold filaments are formed. From an observational point of view the most 
striking consequences are the large oscillations of Lx coupled to the oscillations of the gas density. 
These oscillations may have amplitudes ranging over two orders of magnitude while the X-ray 
surface brightness does not stray dramatically from spherical symmetry. Thus, at least for WOI 
models, the X-ray luminosity does not represent a good diagnostic of the dynamical state of the 
gas when tidal interactions are present. 

The mechanism described above has been proven to be self-sustaining and independent of the 
external cause which initially originated it. In fact, we ran the model MP where, as anticipated 
in §3, the perturber is absent, the galaxy is at rest in a inertial reference system, and a random 
large-scale "noise" is superimposed on the initial velocity field. In model MP the maximum 
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allowed relative amplitude of the perturbation with respect to the unperturbed velocity field is 20 
per cent. Apart from this, the initial conditions are the same as for model M3. The mass budget 
and Lx evolution of model MP are shown in Fig. 8, which should be compared with Fig. 3cf. As 
can be seen, this model behaves mostly as model M3. The spikes in Lx are somewhat larger, but 
this is due to our (arbitrary) choice of the maximum amplitude of the velocity perturbations. We 
point out that model MP is rather similar to models discussed by Kritsuk et al. (1998). However, 
a significant comparison between these models and model MP cannot be made because the latter 
authors stop their simulations after only 50 Myr, a time interval too short to predict the behaviour 
of the flow on cosmological times. The major insight derived from model MP is that the evolution 
of perturbed WOI models is weakly dependent on the specific nature of the perturbation. 

We also ran a pure cooling flow model with a tidal perturbation (model MC): in this case the 
Lx oscillations are significantly reduced (approximately of a factor of ten) with respect to those 
exhibited by model M3. The most important consequence is that tidal interactions cannot be at 
the origin of the spread in the Lx — L-q diagram if the pure cooling flow scenario applies. 

As discussed in §2, in our simulations we neglected the cluster tidal field, which is present in 
equation (5). Although the magnitude of the cluster tidal field during the encounter is negligible 
(see §2.1.1), one can ask whether the cumulative effects of such field over cosmological times are 
significative. In order to check this possibility, as anticipated in §3 we ran the model MT, in which 
-^SN = 0, no perturber is present, and only the cluster tidal field is considered. The resulting 
evolution of this model (not shown here) is very similar to that of model MC, but the amplitudes 
in the Lx oscillations are much reduced. We thus conclude that the cluster tidal field can be safely 
neglected when studying the effects of galaxy-galaxy encounters like those considered here on the 
X-ray properties of early-type galaxies. 

Moreover, it must be added that the role of the cluster tidal field is further reduced by the 
fact that during a Hubble time a galaxy may experience several encounters of the type discussed 
here. A simple estimate of this number can be obtained as follows. Suppose that the galaxy 
belongs to a typical rich cluster like that described in §2.1.1, and containing N'y = 300 galaxies, 
whose luminosity function is given by 

<P{L)dL = N, (^\ exp {-L/L,)d (^) , (15) 

where L^ ~ 4.74 x 10"'^'' Lq, a = 5/4 (Schechter 1976, Sarazin 1986), and with lower and upper 
limits of 10^ -Lq and oo, respectively. Using simple geometrical arguments it is easy to show that, 
for an impact parameter less than or equal to r™"\ the mean number of interactions per galaxy 
during an Hubble time in is given by 

In the equation above k is the fraction of the total number of galaxies with luminosity greater than 
Lb- For example, k ~ 0.01 for Lb = 5 x 10^^ Lq (the luminosity of our galaxy model); k ~ 0.1 for 
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Lb = 1.6 X W^^Lq (one third of the luminosity of our galaxy model). It turns out that ni ~ 2 
for K = 0.01. The same result is obtained considering galaxies with Lb > 1.6 x IO^'^Lq in a poor 
cluster (or loose groups), with Nj' = 30, r^ = 200 kpc and a^i = 250 km s^^ (see, e.g., Bahcall 
1998). 

Obviously the above analysis does not apply to small compact galaxy groups containing only 
a few members. In fact, these systems are highly collisional, and the mean separation between 
galaxies is comparable to their dimensions (see, e.g., Kelm, Focardi, & Palumbo 1998). As a 
consequence, strong gravitational interactions are frequent, the evolution of gas flows is higly 
perturbed, and the results of our simulations can not be directly applied in this case. 

From a more observational point of view, the importance of the discussion above is also 
strengthen by the fact that, while it is true that the majority of elliptical galaxies are found in 
clusters, the galaxies from which the Lx — Lb diagrams are constructed are either in groups or 
in the Virgo cluster (which is less rich than the cluster assumed in our simulations). Hovewer, 
from the discussion above it results that the "turbulent" mechanism here described should be 
quite common in ellipticals in clusters and groups, while galaxies in compact groups can hardly 
avoid strong gravitational interactions, that should produce either strong variations in their X-ray 
luminosity or even substantial degassing. 

As pointed out in §2.2, with the computed models available, we are now in the position to 
evaluate the error induced by neglecting the non-inertial forces associated with the rotation of 
the adopted reference frame. From an analysis of the model's velocity field at different times, it 
turns out that, when the perturber is closest, the Coriolis and centrifugal forces are larger than 
the gravitational forces only within a small volume (with a characteristic size of ~ 2 kpc) at the 
edge of the galaxy facing the perturber, where the total gravity is very low. Later (t — tm ~ 200 
Myr), the gravity largely prevails (by a factor 100-1000) all over the computational grid. This is 
mainly due to the fact that for the assumed orbit ||!7|| oc {t — tm)^^ (see equation [10]), and so it 
quickly reduces after the encounter; moreover, the non inertial force associated with f2 vanishes 
when Tp = r™'", and, for |t — tm\ > r™"^/!/ ~ 100 Myr, \\f2\\ oc |t — tm|~^- As a consequence, we 
believe that our results are reliable, and, although neglecting the fictitious forces due to rotation 
may prevent from obtaining the exact values of the Lx spikes, the essential physics is captured by 
the computed models. 

In conclusion, the main results of this work can be summarized as follows: 

• In WOI models the Lx oscillations induced by tidal encounters last for several Gyrs after the 
encounter is over. The qualitative evolution of these oscillations is not strongly dependent on 
the flow phase (wind, outflow, or inflow) occurring when the encounter starts; moreover, they 
appear to be rather independent of the physical process generating them. The characteristic 
period of such oscillations is of the order of the sound crossing time of the hot gas in the 
galactic inner region. If the encounter happens quite early, when Lsn is still relatively 
high, then a catastrophic drop of Lx is produced. These oscillations can contribute to the 
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observed spread in the Lx — L^ diagram for galaxies in clusters and groups. Of course, this 
mechanism cannot be at the origin of the Lx scatter in isolated galaxies. 

• When compared to WOI models, the amplitude of the Lx oscillations in cooling flow models 
(with LsN = 0) is substantially reduced. As a consequence, tidal interactions cannot be at 
the origin of the observed spread in the Lx — Lb diagram in the "pure" (no SNIa heating) 
cooling flow scenario. 

• In all computed models, the X-ray surface brightness Ex does not appear to be strongly 
distorted by tidal interactions. However, the (angular mean) radial profile of Sx is steeper 
in the inner region in coincidence with the Lx spikes. 

• We found that the tidal stripping due to the perturber is effective. In fact, approximately 
half of the gas content of the galaxy is lost in perturbed WOI models when compared to the 
unperturbed one. In our pure cooling flow model the fraction of lost mass is reduced to one 
fourth. 

• The cluster tidal field alone appears to be ineffective in producing large Lx oscillations when 
compared to tidal interactions. Moreover, during a Hubble time, galaxies in clusters suffers 
several encounters of the kind described here, and so the importance of the cluster tidal field 
on the X-ray properties of early-type galaxies is further reduced. 
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A. The Hydrodynamical Equations with General Source Terms in Non— Inertial 

Reference Systems 

In Appendix A.l we present the treatment of mass, momentum, and energy source terms in 
hydrodynamics, in the general anisotropic case. In Appendix A. 2 the hydrodynamical equations 
with general source terms are expressed in a non inertial reference frame. 



A.l. The Source Terms 

In the following treatment we assume the presence in the (inertial) reference system S* of a 
single source field, the generalization to the case of more than one field being straightforward. 
X = XjCj and v = WjCj are the position and the velocity vectors in S, respectively. Let the 
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source field be described by a distribution function / = /(x, v;t) (see, e.g., BT), so that 
n(x; t) = /jj3 f d^v represents the number density of sources. 

The mass return associated with the source field is given by the function m = m{x, v , n; t) , 
where n = njej is a unitary vector accounting for the possibihty of anisotropic mass sources; the 
physical units of m are mass for unit time for unit solid angle. The total mass return per unit time, 
volume and solid angle associated with the source field is then given by //(x, n;f) = /^jsm/d^v. 
Finally, integrating over the whole solid angle, we obtain the total mass return rate per unit time 
and volume at x: 

M{x;t) = I fid^n. (Al) 

J47r 

For example, for mass sources independent of the source velocity v, /u(x, n; t) = nm; moreover for 
an isotropic source field, A4(x; t) = Aimm. 

The momentum return associated with the source field is given by the vectorial function 
p = m(x, V, n;t)[v + iis(x, v, n; t)n], where u^ is the modulus of the velocity of the material 
associated with the source field along n with respect to the source velocity v. The total 
momentum per unit time, volume and solid angle associated with the source field is then given by 
7r(x, n; t) = j'^-j p / d^v. Finally, integrating over the whole solid angle, 

P(x;t) = / TTci^n. (A2) 

In particular, for mass sources and ejection velocity independent of the source velocity v, 
7r(x, n; t) = nm^ + Ugii], where 

n{x;t)V{yL;t) = f v f d^v (A3) 

is the source streaming velocity field, v = zJjej (see, e.g., BT). Moreover, if m and Ug are also 
isotropic, P(x;t) = A^ v. 

The energy source associated with the source field is made of the contributions of two 
distinct parts, i.e., the internal energy source and the kinetic energy source. The internal energy 
source is described (in analogy with the mass source) by the function e = e(x, v, n;t). The 
physical units of e are energy per unit time per unit solid angle per unit mass. The total internal 
energy per unit time, volume and solid angle associated with the source field is then given by 
e(x, n; t) = /s^3 me f d^v. Finally, integrating over the whole solid angle. 



£{x;t) = I ed^n. (A4) 

J47r 



In particular, for mass and internal energy sources independent of the velocity source v, 

e(x, n; t) = /ue. Moreover, in case of isotropy, <S(x; t) = Me. The kinetic energy source is given by 

k = ir7T,(x, V, n; t)||v + Us(x, v, n; t)n|p. The total kinetic energy per unit time, volume and solid 
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angle associated with the source field is then given by «;(x, n;t) = /^js k f d^v. Finally, integrating 
over the whole solid angle, 

/C(x;t) = / nd'^n. (A5) 

In particular, for mass sources and ejection velocity independent of the velocity source v, 
k(x, n;t) = nr7T,[||v||^ + Tr(cj^) + Ug + 2us < n, v >]/2, where Tr((T^) is the trace of the velocity 
dispersion tensor associated with the source motion: 

«(x; t)afj{x; t) = (vi - Vi) {vj - Vj) f d^v. (A6) 

Moreover, if m and Us are also isotropic, /C(x; t) = A^[||v|p + u^ + Tr(cr^)]/2. Note that a source 
field is called isotropic only if the functions m, p, u^ and e are all independent of n. 



A. 2. The Hydrodynamical Equations in a Non Inertial Reference System 

The basic hydrodynamical equations with general source terms are here derived (in the 
inviscid case) by using the standard approach based on the Reynolds Transport Theorem. In fact, 
this formulation is particularly useful when moving to a non-inertial reference frame (see, e.g., 
Narasimhan 1993). 

As usual, we denote with D/Dt = d/dt + Vid/dxi the Lagrangian derivative when expressed 
in Eulerian form. After some manipulation of the conservation laws expressed in integral form, 
and application of the Transport Theorem, one obtains the continuity equation 

^ + pVx-u = A^. (A7) 

The three components of the momentum equation are given by 

p-^ = pg-VxP + P-A^u, (A8) 

where g is the total physical acceleration and p = p(x; t) is the thermodynamical pressure. Finally, 
the energy equation is derived: 

^ + (^ +p)Vx • u = £: + /C + ^||u||2- < P,u > -C. (A9) 

where E is the internal energy per unit volume, and C describes the radiative losses per unit 
volume and unit time. Note that thermal conduction is not considered here. 

In case of a completely isotropic source field (the case of our model galaxies), the momentum 
and energy equations can be rewritten as: 

p-^ = pg-VxP + A^(v-u). (AlO) 
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DE 



+ {E+p)V^-u 



M^ 



u 



■v||2+7W 



2 2 



C. 



(All) 



We move now from S* to a non-inertial frame S'. As well known, the position, velocity, and 
acceleration vectors x, x, x (in S) are related to the corresponding vectors in S' by 



X = R + Ol, 



X 



X 



A + 0[^ + 2f2A^+f2A^+f2A{f2A$,)], 



(A12) 

(A13) 

(A14) 



where R(t), V(t) and A(t) are the position, velocity, and acceleration of the origin of S" (see, e.g., 
Arnol'd 1980). 0{t) G SO (3) is a rotation matrix, and f2 is the dual of the antisymmetric matrix 
O^O. Geometrically, fi is the angular velocity of S' with respect to S resolved along the basis of 
S' . When resolved in S", we have u) = Ofl. 

Now, the scalar continuity and energy equations are transformed in the corresponding 
equations in S' . First, the relation between the velocity field u(x;t) (in S) and i'{$,]t) (in S') is 
obtained from the definition of velocity field as Lagrangian derivative of the position vector, and 
so from equation (A13) u(x;t) = V + 0[v{$,;t) + fl A ^. Second, in S", the expression for the 
Lagrangian derivative is D/Dt = d/dt + i^id/d^i, moreover, it is easily proved that Vx • u = Ve • v. 

The transformation of continuity and energy equations is now immediate; some care is 
required for the momentum equation, when considering the Lagrangian derivative of the velocity. 
The final equations (for isotropic sources) in conservative form are 



^+V^-pu = M, 



dpu 



+ V^- (pz/0iv) = pO^{g- A) - p[2QAu+ flA£, + flA (^A^)] -V^p + Mw\ 



dE 
'dt 



M, 



+ V^- Ei'= -pV^-u+ — 1|«^- V H +M 



ul Tr(cj2) 

e + — H ^ — - 

2 2 



C. 



(A15) 
(A16) 

(A17) 
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Fig. 1. — The modulus of the dimensionless radial (solid line) and tangential (dashed line) 
component of the cluster tidal field, as a function of the cluster radius. The normalization constant 
is |Tp|™™ = GMp/rp, for Vp = 100 kpc. The explicit expression of the q function entering equation 
(8) can be found in Ciotti & Dutta (1994). The cluster mass is obtained by using the standard 
relation 47rG/9ci(0)r^i = 9a'^^. 

Fig. 2. — Logarithmic density distribution and velocity field of the gas in model Ml at different 
times. The perturber is on the right. 

Fig. 3. — The upper panels show the evolution of Lx for the models Ml, M2, and M3 (solid lines), 
for the model MO (dashed lines), and Lsn (dotted lines). The two dots refer to the times indicated 
in Fig. 2bc, and Figs. 6 and 7. The evolution of the total gas mass (heavy lines), hot gas mass 
(medium lines), and cold gas mass (light lines) are shown in the lower panels. Dashed lines refer 
to model MO, sohd lines to models Ml, M2, and M3. 

Fig. 4. — Logarithmic density distribution and velocity field of the gas in model Ml near the 
galactic center at two different times (in particular the left panel is the enlargment of the central 
region of Fig. 2b). The turbulent regime of the flow, and the transient cold filaments are apparent. 

Fig. 5. — Time evolution of Lx for the cooling flow, tidally perturbed MC model (solid line). Lx 
of the same model without any perturbation (dashed line) is shown for comparison. 

Fig. 6. — X-ray surface brightness Sx of model Ml at t = 9.9 Gyr (solid curves) and t = 11.1 Gyr 
(dotted curves). These times are the same as in panels b and c of Fig. 2. The perturber is located 
on the right. Labels are in kpc. 

Fig. 7. — Radial distribution of the angular averaged surface brightness Sx shown in Fig. 6. Solid 
and dotted lines correspond to a Lx maximum and minimum, respectively (see the black dots in 
Fig. 3a). Te = y^rvTd — 4.9 kpc is the effective radius of the stellar density distribution. T,^^^ is the 
central value of the X-ray surface brightness of model MO. 

Fig. 8. — X-ray luminosity Lx and gas mass evolution for the model MP. In the upper panel the 
dotted line represents LgN . In the lower panel the evolution of the total gas mass (heavy line) , hot 
gas mass (medium line), and cold gas mass (light line) is shown. 
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